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ABSTRACT 

We present a kinematical study of 29 spiral galaxies included in the Spitzer Survey of 
Stellar Structure in Galaxies, using Hct Fabry-Perot data obtained with the Galaxy 
Ha Fabry-Perot System instrument at the William Herschel Telescope in La Palma, 
complemented with images in the R band and in Ha. The primary goal is to study 
the evolution and properties of the main structural components of galaxies through 
the kinematical analysis of the FP data, complemented with studies of morphology, 
star formation and mass distribution. In this paper we describe how the FP data 
have been obtained, processed and analysed. We present the resulting moment maps, 
rotation curves, velocity model maps and residual maps. Images are available in FITS 
format through the NASA/IPAC Extragalactic Database and the Centre de Donnees 
Stellaires. With these data products we study the non-circular motions, in particular 
those found along the bars and spiral arms. The data indicate that the amplitude of 
the non-circular motions created by the bar does not correlate with the bar strength 
indicators. The amplitude of those non-circular motions in the spiral arms does not 
correlate with either arm class or star formation rate along the spiral arms. This 
implies that the presence and the magnitude of the streaming motions in the arms is 
a local phenomenon. 
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1 INTRODUCTION 


Galaxies are the basic building blocks of the Universe and 
understanding their nature, formation and evolution is cru- 
* Email: serroz@iac.es cial to many areas of current astrophysical research. Con- 
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cretely, nearby galaxies are the vivid result of the evolution 
of the Universe, they contain the footprints of the evolution 
processes that have led to the present status of the galaxies. 
To better understand the evolutionary history of the Uni¬ 
verse, it is necessary to study the structure and dynamics of 
current galaxies. By studying the kinematics of galaxies we 
can trace the motions of both their baryonic (gas and stars) 
and dark matter. 

Radial velocity measurements have been traditionally 
performed with radio observations, mainly using the 21-cm 
Hi line. Tracing this neutral hydrogen line lets us study most 
of the gas content of the galaxy, usually traced out to three 


or four times beyond the visible disc. 

Rots 

(19751; 

Bosma 

(1978 

); van der Hulst (1979); Bosma 

(1981); Gottesman 


(1982) and others demonstrated the power of rotation curves 


in deriving the total mass distribution of disc galaxies. The 
drawback with these Hi radio observations was the poor an¬ 
gular resolution. The highest angular resolution (~ Q") 21- 
cm Hi surveys of nearby galaxies to date have been carried 
out by the Hi Nearby Galaxy Survey (THINGS; [Walter et al. 


2008), using the Very Large Array (VLA) of the NRAO. GO 


radio observations have also been used traditionally to trace 
the molecular gas. Nowadays, GO observations provide high 
angular resolution (comparable to optical), and also high 
spectral resolution (from one to several km s“^), but in very 
small fields. 

Optical (Ha, Nil) observations yield high arcsecond an¬ 
gular resolution data. Ha is often the brightest emission line 
in the visible wavelength range due to the cosmic abundance 
of hydrogen. In spiral galaxies, this line traces primarily the 
ionized gas in Hll regions around young massive stars. In the 
20th century, most of the optical observations were based 
on long-slit spectroscopy. Long-slit observations have been 
traditionally used to deduce the rotation curves of galaxies 


Rubin et al.|I985| 

[Mathewson, Eord & Buchhorn|I992 

sic & Salucci|I995| 

[Mathewson & Eord|I996||Gourteau 


jPer^j 




However, 3D spectroscopy [IFU, Fabry-Perot (FP), 
multi long-slit spectroscopy] of galaxies is one of the best 
methods to obtain detailed information on the kinematics of 
galaxies. Velocity maps derived using 3D spectroscopy repro¬ 
duce the complete velocity field, contrary to longslit spectra. 
Ha observations using FP instruments have been used for 


(Tully 1974 Deharveng & Pellet 

1975 

or 

de Vaucouleurs & Pence||I980 

); and 


have been designed to derive velocity fields of nearby galax¬ 
ies (e.g., Atherton et al.|I9^ or jBoulesteix et al.|I984 ), and 


to create high signal-to-noise ratio (S/N) rotation curves for 


many spiral galaxies (e.g.. 

Marcelin, Boulesteix & Georgelin 

I985||Bonnarel et al.|I988 

[Pence, Taylor & Atherton|I990[ 

Gorradi & Gapaccioli||I99; 

[ [Amram et al.[[I992[ [Gecil, Wil- 

son & Tully||I992| [Amram et al.[[I994[ [Sicotte, Garignan & 

Durand|I996[ [Ryder et al.|I998[[ Jimenez-Vicente et al.[I999[ 

Knapen et al.[[2000[and m^ 

my more). 

presented an Ha EP survey of 

Ghemin et al.[ (2006) 


an Ha FP survey complementary to the Spitzer Infrared 
Nearby Galaxies Survey (SINGS; [Kennicutt et al.|[2003l l. 
One of the largest and most important FP surveys in Ha 
is the Gassendi HAlpha survey of SPirals (GHASP; Garrido 
et al.|20Q'2 ). The GHASP survey consists of a sample of 203 


spiral and irregular galaxies that have been observed with 
a sampling about 16 km s~^ in velocity and an average 3 


arcsec in angular resolution (see Epinat, Amram V Marcelin 


|2008| GHASP VH hereafter, for a complete list of data and 
resolutions). Ha FP spectrographs are nowadays being used 
for studying the kinematics of several kinds of galaxies, e.g., 
bulgeless galaxies ( Neumayer et al.]|20II ), star burst galax¬ 
ies (Blasco-Herrera et al. 2013), or interacting galaxies in 
compact dwarfs ( Torres-Flores et al.||20I4 ). 

One of the newest FP spectrographs is the Galaxy Ha 
Fabry-Perot System (GHaFaS). It is a visiting instrument 
on the William Herschel Telescope (WHT) in La Palma. The 


GHaFaS instrument has been operative since 2007 (Fathi 
et al.|[20Q7 ), and has been used to study pattern speeds of 


bars and spiral arms ( Fathi et al.|2009 ), to analyse the kine¬ 
matics of planeta ry nebulae (| Santander-Garcfa et al.|20I0 ), 
and the gas flows (Font et al.|20IIb), star formation and the 


kinematics of interacting galaxies (Zaragoza-Gardiel et al. 


2013 


2014) and starburst galaxies (Blasco-Herrera et al. 


2013). It has also been used to study the resonance radii 


and interlocking resonance patterns in galaxy discs (Font 
et al.|2011al|2014a[ | 


We have designed an observing programme to obtain 
FP data of 29 spiral galaxies with the GHaFaS instrument, 
as part of the ancillary data of the Spitzer Survey of Stel¬ 
lar Structure in Galaxies (S^G; Sheth et al. 2010). The S^G 
survey has obtained 3.6 and 4.5 /xm images of 2352 nearby 
galaxies using the Infrared Array Gamera (IRAG; [Fazio et ST 


2004). The sample is composed by galaxies that fulfil these 


requirements: d < 40 Mpc, rrih < 15.5, 1)25 > I arcmin, 
and includes galaxies selected using values from HyperLeda 
(Paturel et al. 2003), with a radial velocity Vradio < 3000 
km s“^ and galactic latitude \b\ > 30°. The cornerstone of 
the S^G survey is the quantitative analysis of photometric 
parameters, enabling a variety of studies on secular evolu¬ 
tion, outer disc and halo formation, galaxy morphology, etc. 
The data have been made public^ The first papers resulting 


from S G survey have been summarised injHolwerda et al. 

dMIf . 

The combination of the S^G mid-IR images with these 
complementary kinematic FP data at high resolution will al¬ 
low us to tackle the scientific goals of this study, summarized 
as followings: 


(i) Perform a detailed study of the kinematical interplay 
between the interstellar medium and regions of star forma¬ 
tion, dust, or other activity. 

(ii) Use the kinematics as a probe of the secular evolution, 
as manifested in the structural components such as bars, 
spiral arms, rings, lenses, etc. 

(hi) Exploit the high angular resolution of the kinematic 
data to study the inner parts of rotation curves, and relate 
the galaxy kinematics to the mass distribution and specific 
observed properties, as well as probing the coupling between 
the stellar density and the gravitational potential in the in¬ 
ner parts of the galaxies. 

(iv) Study possible deviations of the rotation curve that 
are caused by lopsidedness or asymmetries in the outer parts 
of the disc, which would help to pin down the interplay be¬ 
tween dark matter (DM) and stars. 


^ (http://irsa.ipac.caltech.edu/data/SPITZER/S4G/) 
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(v) Explore the outer disc kinematics, to establish 
whether discs are generally cold and thus to constrain the 
halo properties. 


This data paper is the second of a series, and presents 
the complete set of FP and narrow-band imaging data for 


the 29 galaxies in our kinematical study. In Paper I (Erroz- 


Ferrer et al. 2012), we illustrated the data and methods. 


and discussed the kind of results that the main survey will 
provide through a detailed analysis of NGC 864, an archety¬ 
pal barred spiral galaxy. In Paper III (Erroz-Ferrer et ah; in 
prep.) we will present the results from a study of the inner 
part of the rotation curves of the galaxies of the sample; 
and Paper IV (Leaman et ah; in prep) will study the outer¬ 
most parts of the discs of the galaxies of the sample and the 
relationship between their kinematics and DM. 

This paper is organized as follows: Section ^ gives a 
description of the sample selection, and Section |3 describes 
the observations. Section |4] describes the data reduction and 
results. The derived rotation curves and non-circular mo¬ 
tions maps are presented in Sections [5.1|and | 5.2| respectively. 
These results are discussed in Section |5.4| and Section 
presents our conclusions. 


2 TARGET SELECTION 


The galaxies in our survey satisfy the following require¬ 
ments: first of all, the declination should be higher than 
-10° so that the altitude in the sky could be enough at the 
time of the observation in La Palma, and only galaxies with 
inclinations between 0° and 70° were selected. 

According to the instrument specifications, the galaxy 
should fit well in the GHaFaS FOV of 3.4 x 3.4 arcmin. 
Therefore, galaxies with diameters between 2 and 3.4 ar¬ 
cmin were selected. Secondly, the range in velocity of the 
galaxy should not be higher than the free spectral range 
(FSR) of the instrument (see Sect. for details), although 
we did not have this ancillary information for all the galax¬ 
ies. Furthermore, the data reduction processes require that 
the galaxy should have at least three bright point sources in 
the field, i.e. three stars or two stars and the nucleus. 

The observed sample was selected at the time of the 
observations, looking at the visibility on the sky, but also 
requiring a spread in morphological type, bar presence and 
finally, the availability of ancillary data (interferometric GO, 
Hi, ultraviolet or Spitzer mid-IR). 

The final sample consists of 29 galaxies. For the morpho¬ 
logical classification, we have followed the up-to-date Mid- 
IR Glassifications for S^G Galaxies in the Gomprehensive 
de Vaucouleurs revised Hubble-Sandage System (GVRHS, 
Buta et 'M^|2015 ). The morphological classification also in¬ 
cludes bar presence (9 SA galaxies, 12 SAB galaxies and 7 
SB galaxies). In Fig.[^ we present the morphological classifi¬ 
cation as a histogram. The general properties of the galaxies 
in the sample are presented in Table ^ As noted in the foot¬ 
note of Tablethe galaxy NGG 7241 is nearly edge-on. The 
galaxy has an unusual interacting companion in the line of 
sight, its photometric inclination is biased much lower and 
therefore we do not include it in the analysis part of this 
paper. A detailed analysis of it is presented separately in 
Leaman et al. (in prep.). 


MORPHOLOGICAL TYPE 



II SB 
■ SAB 
= SA 


Figure 1. Morphological distribution of the galaxies of the sam¬ 
ple, organised by morphological type and also by bar presence. 


3 OBSERVATIONS 


The Fabry-Perot observations were carried out using the 
GHaFaS instrument mounted on the 4.2-m WHT in La 
Palma. The observations were done during 24 nights be¬ 
tween September 2010 and March 2013. 

The instrument is a Fabry-Perot interferometer that 
provides high spectral resolution and seeing-limited angu¬ 
lar resolution data, within a 3^4 x 3^4 FOV. The instru¬ 
ment comprises a focal reducer, a filter wheel, a Fabry-Perot 
etalon, an image photon-counting system (IPGS), and a cali¬ 
bration lamp (neon source). The etalon used has an interfer¬ 
ence order of 765 and a FSR of 391.9 kms“^ (8.6 A). Each 
observational cycle consists of 48 steps (named channels) 
through the etalon for 10 seconds each. The number of cy¬ 
cles depends on the galaxy magnitude, but usually amounts 
to some 3 hours (corresponding to 22 cycles). However, due 
to weather and instrumentation problems, not all the galax¬ 
ies were observed for the same length of time (all the infor¬ 
mation about the observations is collected in Table |^. The 
high spectral resolution mode was used, achieving a velocity 
sampling of ~8 km s~^ with a pixel scale of 072 in IKxlK 
pixel images. In Table we summarize the instrumental 
setup used in the observations. 

To perform flux calibration, Ha narrow-band images 
were observed together with the Fabry-Perot observations, 
using the Auxiliary port GAMera (AGAM), permanently 
mounted at a folded-Gassegrain focus of the WHT (Benn, 
Dee V Ag6cs|200'8 ). The FOV in imaging mode is 8 arcmin. 


with a pixel size of 0.25 arcsec. The R-band images were 
taken using a Bessel filter with central wavelength 6500A 
and full width half maximum (FWHM) of 1360A, and with 
a Sloan 6228/1322 filter (central wavelength/FHWM, in A). 
We used four different Ha filters depending on the galaxy’s 
redshift (6577/23, 6589/24, 6613/24 and 6631/17; central 
wavelength/FHWM, in A). The exposure times per galaxy 
were 3x20 seconds for the R-band images and 3x100 sec¬ 
onds for the Ha images. The typical seeing was around 1 
arcsec, ranging from 076 to l76. 
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Table 1. General properties of the galaxies in the sample. Notes. (1) Updated morphological classifications from the CVRHS, where 
“double stage” galaxies are allowed (i.e., large-scale SO or SO/a galaxies with smaller-scale inner spirals) (2) Adopted values of the 
distances, calculated after applying the Virgo, GA and Shapley corrections, with Hq= 73 ± 5 km s“^ Mpc“^, from the NASA/IPAG 
Extragalactic Database (NED). The uncertainties in the distance measurements have been adopted as 20 per cent of the value. (3) B 
magnitude from The Third Reference Gatalogue of Bright Galaxies (RG3; |de Vaucouleurs et al.|p~991| >. (4) Absolute B magnitude 
measured using the distance and m^ of columns III and IV. (5) Absolute 3.6/rm magnitude measured using the distance of column III 
and asymptotic magnitude at 3.6 jim from the ellipse fitting to the 3.6 /rm S'^G images (Munoz-Mateos et al. in prep). NGG 7241 
is nearly edge on, but due to a line of sight companion galaxy, the photometric inclination is biased much lower and should not be 
considered representative. 


Galaxy 

■ name 

mid-IR morphology 
(1) 

d(Mpc) 

(2) 

niB 

(3) 

Mb 

(4) 

Ms.6 
(5) 

NGG 428 

SAB(s)dm 

15.9 ± 3.2 

11.95 

-19.06 

-19.20 

NGG 691 

(R)SA(rs,rl)ab 

35.7 ± 7.1 

12.28 

-20.48 

-21.36 

NGG 864 

SAB(rs)bc 

20.9 ± 4.2 

11.62 

-19.98 

-20.44 

NGG 918 

SAB(s)cd 

20.5 ± 4.1 

13.07 

-18.49 

-20.24 

NGG 

1073 

SB(rs)cd 

I6.I ± 3.2 

11.68 

-19.35 

-19.75 

NGG 

2500 

SAB(s)cd 

9.8 ± 2.0 

12.22 

-17.73 

-18.14 

NGG 

2541 

SA(s)dm 

10.4 ± 2.1 

12.25 

-17.84 

-17.90 

NGG 

2543 

SAB(s,bl)b 

37.4 ± 7.5 

12.94 

-19.92 

-20.96 

NGG 

2712 

(R')SAB(rs,nl)ab 

29.5 ± 5.9 

12.78 

-19.57 

-20.67 

NGG 

2748 

(R')SAB(rs)bc 

25.1 ± 5.0 

12.39 

-19.61 

-20.84 

NGG 

2805 

(R)SA(s)c pec 

28.7 ± 5.7 

11.79 

-20.50 

-20.51 

NGG 

3041 

SA(rs)bc 

23.9 ± 4.8 

12.30 

-19.59 

-20.52 

NGG 

3403 

SA(rs)c: 

22.8 ± 4.6 

12.94 

-18.85 

-19.63 

NGG 

3423 

SA(s)bc 

14.3 ± 2.9 

11.60 

-19.18 

-19.58 

NGG 

3504 

(R'^)SAB(rs,nl)a 

27.8 ± 5.6 

11.62 

-20.60 

-21.65 

NGG 

4151 

SABa(l,nl)0/a 

20.0 ± 4.0 

11.36 

-20.15 

-21.40 

NGG 

4324 

(L)SA(r)0+ 

13.6 ± 2.7 

12.50 

-18.17 

-19.51 

NGG 

4389 

SB(rs)a[d] 

13.8 ± 2.8 

12.55 

-18.15 

-19.11 

NGG 

4498 

SB(rs)d 

I4.I ± 2.8 

12.77 

-17.98 

-18.51 

NGG 

4639 

(R')SAB(rs,bl)ab 

13.9 ± 2.8 

12.19 

-18.53 

-19.39 

NGG 

5112 

SB(s)cd 

20.2 ± 4.0 

12.63 

-18.90 

-19.49 

NGG 

5334 

SB(rs,xir)cd 

24.2 ± 4.8 

12.88 

-19.04 

-20.03 

NGG 

5678 

(R'L)SA(rs)b: pec 

33.3 ± 6.7 

12.02 

-20.59 

-21.87 

NGG 

5740 

(RL)SAB(r)ab 

27.0 ± 5.4 

12.60 

-19.56 

-20.65 

NGG 

5921 

SB(rs)b 

26.2 ± 5.2 

11.68 

-20.41 

-21.31 

NGG 

6070 

SA(rs,nrl)c 

33.6 ± 6.7 

12.42 

-20.21 

-21.70 

NGG 

6207 

SAB(rs)cd 

18.5 ± 3.7 

11.86 

-19.48 

-19.84 

NGG 

6412 

SB(rs)cd 

23.7 ± 4.7 

12.38 

-19.49 

-20.07 

NGG 

7241 

Scd sp / E(d)7 

22.4 ± 4.5 

13.23 

-18.52 

-20.51 


Table 2. Journal of observations. 


Observations 

Telescope 

4.2-m WHT 


Fabry Perot 

Instrument 

GHaFaS 



Spatial sampling 

FOV 

Pixel scale 

202” X 202” 

0'.'2 


Galibration 

Neon comparison light 

A 6598.95 A 


Gharacteristics @ Ha 

Interference order 

FSR 

Finesse 

Spectral resolution 
Instrumental FWHM 

765 

391.9 km s-i (8.66 A) 

20 

18179 

19.6 km s“^ 


Spectral sampling 

Detector 

Number of scanning steps 
Sampling step 

IPGS 

48 

8.2 km s-i (0.18 A) 

Imaging 

Instrument 

AGAM 



Spatial sampling 

Galibration 

FOV 

Pixel scale 

Standard stars 

8 arcmin 

0'/25 
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Table 3. Log of the observation. Notes: (1) The adopted format for the filters is central wavelength/FHWM (2) The seeing was measured 
from the ACAM images. To perform the continuum subtraction, the Ha image or the R-hand image were degraded to the same spatial 
resolution (the worst seeing of the two images), and this number has been adopted as the final seeing. 


Galaxy name 

Date 

FP texp 

(s) 

FP Hq filter(i) 

(A/A) 

ACAM Hq filter(i) 

(A/A) 

seeing^^^ 

(”) 

NGC 428 

Sept. 2010 

3840 

6580/23 

6589/24 

1.0 

NGC 691 

Jan. 2012 

11040 

6623/23 

6631/17 

1.2 

NGC 864 

Sept. 2010 

13920 

6600.5/25 

6589/24 

0.9 

NGC 918 

Sept. 2010 

2880 

6600.5/25 

6589/24 

0.9 

NGC 1073 

Sept. 2010 

20640 

6580/23 

6589/24 

0.9 

NGC 2500 

Jan. 2012 

10560 

6580/23 

6577/23 

1.0 

NGC 2541 

Jan. 2012 

10560 

6580/23 

6577/23 

1.2 

NGC 2543 

Jan. 2012 

12960 

6623/23 

6613/24 

1.3 

NGC 2712 

Feb. 2012 

10560 

6600.5/25 

6613/24 

1.4 

NGC 2748 

Jan. 2012 

9600 

6600.5/25 

6589/24 

1.4 

NGC 2805 

Mar. 2013 

9120 

6597.5/17.6 

6613/24 

1.1 

NGC 3041 

Jan. 2012 

10080 

6600.5/25 

6589/24 

1.4 

NGC 3403 

Jan. 2012 

10560 

6580/23 

6589/24 

1.1 

NGC 3423 

Feb. 2012 

9120 

6580/23 

6589/24 

1.5 

NGC 3504 

Jan. 2012 

15360 

6600.5/25 

6589/24 

1.4 

NGC 4151 

Feb. 2012 

13920 

6580/23 

6589/24 

1.2 

NGC 4324 

Mar. 2013 

8160 

6597.5/17.6 

6613/24 

1.4 

NGC 4389 

Jan. 2012 

12000 

6580/23 

6577/23 

0.8 

NGC 4498 

Jun. 2011 

10080 

6598/18 

6589/24 

0.9 

NGC 4639 

Jun. 2011 

8640 

6583/15.5 

6589/24 

1.0 

NGC 5112 

Jun. 2011 

9120 

6583/15.5 

6589/24 

0.6 

NGC 5334 

Jun. 2011 

11520 

6598/18 

6589/24 

1.2 

NGC 5678 

Jan. 2012 

10560 

6600.5/25 

6613/24 

1.2 

NGC 5740 

Jun. 2011 

11040 

6598/18 

6589/24 

1.5 

NGC 5921 

Mar. 2013 

7200 

6597.5/17.6 

6589/24 

1.3 

NGC 6070 

Jun. 2011 

9600 

6608/16 

6613/24 

1.5 

NGC 6207 

Jun. 2011 

8640 

6583/15.5 

6589/24 

1.3 

NGC 6412 

Jun. 2011 

13920 

6598/18 

6589/24 

1.3 

NGC 7241 

Jun. 2011 

7680 

6598/18 

6589/24 

1.1 


4 DATA REDUCTION 
4.1 Ha imaging: ACAM 


The ACAM images were reduced initially using IRAF. First, 
bias and flat corrections were made. Then, the sky was sub¬ 
tracted before and after the combination of the exposures. 
The images were astrometrically calibrated using KOORDS in 
KAPPA first, setting a preliminar astrometry using the DSS 
images as a reference. Afterwards, we used GAIA in STAR- 
LINK to improve the astrometry with the GSC-2 catalogue 
at ESO. The next step was aligning the i?-band and Ha im¬ 
ages for each galaxy, subtracting the continuum using the 


procedures outlined in Knapen et al. (2004) and Sanchez- 


Gallego et al. (2012). Finally, the continuum-subtracted Ha 


images were flux calibrated using spectrophotometric stan¬ 
dard stars observed at the time of the observations. 

The resulting continuum-subtracted image of 
NGC 2748 is presented in panel f) of Fig. The im¬ 
ages of the other galaxies in the sample are presented in 
Appendix A. 


4.2 Fabry-Perot: GHaFaS 

The basic custom reduction of GHaFaS data cubes has been 
explained in previous papers in the literature (e.g. 


Hernan- 


dez et al^|2008 or Blasco-Herrera et al.||2010| ), and 


was in¬ 


troduced in Paper I. The data of all the galaxies have been 
reduced following these steps: de-rotation, phase-correction, 
combination and wavelength calibration, astrometry correc¬ 
tion, spatial smoothing, continuum subtraction, removing 
sky lines, flux calibration and creation of the moment maps. 


4-2.1 De-rotation, phase-correction, combination and 
wavelength calibration 

There is no suitable de-rotator at the Nasmyth focus of the 
WHT. Therefore, the first step is to de-rotate the data cubes. 
To do this, we followed Blasco-Herrera et al. ( 2010| . At least 
two point sources in all the planes are selected (ideally stars), 
and selected throughout the cubes. After that, each plane 
of the subsequent cubes is de-rotated to the same position. 

Due to the nature of the FP data, any 2D spatial trans¬ 
formation (rotation or translation) has to be done also in 
the third dimension (the spectral dimension), and therefore 
de-rotation must be done simultaneously with the phase cali¬ 
bration. The wavelength calibration is also performed at this 
point. Lamp exposures were taken between galaxy observa¬ 


tions and are used to calibrate in wavelength (see Garignan 
et al.|2Q08 ). The integration of all the cubes from the various 
cycles is performed after the de-rotation and phase calibra¬ 
tion. The result is one 48-channel cube, phase-corrected and 
calibrated in wavelength. 
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4 . 2.2 Astrometry calibration 

The data are placed on an astrometrically correct spatial 
grid by comparing the positions of stars in the ACAM Ha 
images. The astrometric calibration of the ACAM images is 
explained in Sect. |4.l[ 


4-2.3 Spatial smoothing 

The instrument delivers data cubes with a spatial scale of 
O'.'2/pix, and therefore the spatial resolution is limited by the 
seeing. If the de-rotation process is not correctly performed, 
the spatial resolution might change (i.e., the addition of mis¬ 
aligned cubes would result in a blurred cube, the PSF of the 
point source rises and therefore the seeing increases). Thus, 
the de-rotation output has been carefully checked in order 
to minimise changes in the resolution. 

We applied a 2D Gaussian smoothing kernel with a 
FWHM of 3 pixels to the data in order to improve the signal- 
to-noise ratio (S/N) without degrading the angular resolu¬ 
tion too much. To apply this Gaussian kernel, we have used 
the IDL task filter_image.pro. 


have dealt with the channels as if they were images to cor¬ 
rect for flat fielding. With flat of GIPSY, we fitted a polyno¬ 
mial to the background, eliminated the gradient across the 
affected planes and flattened the background. 


4-2.6 Flux calibration 


We calibrate the Fabry-Perot data by using calibrated 
narrow-band Ha imaging. As explained in Sect.[^ the kine¬ 
matic data have been observed together with AGAM data. 
The procedure was explained in Paper I: fluxes from se¬ 
lected Hii regions in both AGAM and Ha FP data are com¬ 
pared, and fitted to a linear relationship. As a consequence, 
the zeroth moment maps (intensity maps, see the follow¬ 
ing section) are flux-calibrated. GHASP VH uses a different 
method to flux-calibrate their FP observations. Instead of 
carrying out their own observations, they compare the inte¬ 
grated Ha flux of their galaxies to that of the narrow-band 
observations in James et al. (2004). Although our method 


requires more observing time at the WHT, it is more ac¬ 
curate as both FP and narrow-band imaging observations 
are carried out under the same atmospheric conditions and 
airmass. 


4-2.4 Continuum subtraction 


The resulting de-rotated, phase-corrected and wavelength- 
calibrated data cube is then imported into GIPSY (Groningen 
Image Processing System; van der Hulst et al.||199^ . The 


continuum was subtracted using GONREM of GIPSY, which 
estimates the continuum level from line-free frames. To de¬ 
termine the channels free of line emission, we visually check 
the images and also the spectra, after additional smoothing. 

For NGG 2748, NGG 5678 and NGG 6070, the velocity 
amplitude is equal or larger than the FSR and there are no 
channels free of line emission. The Ha line falls into the next 
interference order and reappears in the first or last channels 
(what we call “peak intrusion”). Those channels containing 
emission from a peak intrusion are copied after or before the 
first or last channel, where they should really be, increasing 
the number of channels per cube. Then, we divided the cube 
spatially in two parts (roughly corresponding to the areas 
where emission from the receding and approaching halves 
is seen). The channels free of emission were identified sepa¬ 
rately for each sub-cube, and the continuum was calculated 
and subtracted separately in each sub-cube. After the con¬ 
tinuum subtraction, the two cutouts have been combined 
again. 


4-2.5 Sky line removal 

Stationary secondary peaks sometimes appear in the data 
cubes, usually located in three channels, peaking in the mid¬ 
dle one and with a circular light gradient along the image 
(more light in the centre and gradually less light as leav¬ 
ing the centre). These are OH sky emission lines due to 
airglow, peaking near the following velocities/wavelengths: 
690 km s“^/6577.3 A (present in NGG 2500, NGG 2541, 
NGG 4389 and NGG 6207), 1080 km s”^ / 6586.5 
A (NGG 428, NGG 1073, NGG 3403, NGG 4151, NGG 4639 
and NGG 5112) and 1480 km s“V6596.7 A (NGG 4498, 
NGG 5334 and NGG 6412). To remove this contribution, we 


4-2.7 Moment maps 


To create the moment maps we have used the MOMENTS 
task of GIPSY, which performs an intensity weighted mean 
of the physical coordinates along the profil^ We imposed 
the condition that the line emission had to be present in at 
least three adjacent channels and at a level above a certain 
noise level a in the prohles. To do that, we determined a for 
each galaxy using STAT in GIPSY, and we created the sets of 
moment maps with emission above 5cr. Also, we imposed the 
condition that secondary peaks are not taken into account 
(although secondary peaks are assumed not to be present 
after the previous reduction processes). 

Subsequently, from the smoothed and continuum- 
subtracted cube we computed the moment maps for each 
galaxy, specifically the moment maps of order zero (inten¬ 
sity map), order one (velocity map along the line of sight) 
and order two (velocity dispersion maps). The velocity dis¬ 
persion maps will be used in Paper HI, and they will be 
corrected for instrumental, thermal and natural li ne broad- 
ening. The natural width has a value of 3 km s“^ (jO’dell & 


Townsle;^|1988 ). The thermal width corresponds to 9.1 km 
s“^, assuming a temperature of 10^ K (Osterbrock & Fer- 
land|2Q06 ). The instrumental width for each galaxy was ob¬ 
tained from the data cube taken with the calibration lamp, 
following the procedures in ?, and is 8.3 km s“^. 

After creating the moment maps, the stars and central 
regions may leave a residual (i.e., the continuum may not be 
properly determined and removed), and therefore the out¬ 
puts from MOMENTS need to be masked. Specifically, for the 
AGN hosts NGG 4151 and NGG 4639, the central emission 
saturated the detector and in consequence the continuum 
subtraction could not be removed properly, leaving a resid¬ 
ual that had to be masked out. Also, we have masked the 
central regions in NGG 4324, NGG 5740 and NGG 5921, 


^ https://www.astro.rug.nl/~gipsy/tsk/moments.del 
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since the signal there is due to continuum emission (i.e., no 
line emission is detected). 

For each galaxy, we created a six-plots figure, present¬ 
ing the resulting velocity maps and the outputs from the 
analysis. In Fig. we present the resulting moment maps 
for NGC 2748 in the top images. The hgures for the other 
galaxies in the sample are collected in Appendix A. 


5 DATA ANALYSIS 
5.1 Rotation curves 

To understand the circular motions within a galaxy, rota¬ 
tion curves have been used widely. Here, we study the ro¬ 
tation of the gas, and in particular, the gas from the Hll 
regions that have been heated by massive, young stars. To 
extract the rotation curve from the velocity map we have 
used the ROTCUR task in GIPSY. This procedure is based on 
the tilted-ring method described by |Begeman (1989), where 
each ring can be defined by the parameters: inclination (i), 
position angle (PA), centre (xq and yo) and systemic veloc¬ 
ity (I’sys)- If we assume that the radial and vertical veloc¬ 
ities are negligible, the observed velocity can be expressed 
as Vobs{R,0,i) = Vsys + Vrot{R,0) cos 0 sin i, where 0 is the 
azimuthal angle in the plane of the galaxy and depends on 
i, PA, xo and yo. 

It is not easy to determine the rotation curve parame¬ 
ters, as deviations from rotational motions can be present 
inside the galaxies, some of our galaxies are not much in¬ 
clined, and most of the velocity fields are sparsely sampled. 
However, we have proceeded to ensure that we are obtain¬ 
ing the true rotation curve while minimising possible errors. 
In order to determine the starting values (i.e., initial condi¬ 
tions) of z, PA, Xq, yo and I’sys, our first step was to search 
the literature (RC3, de Vaucouleurs et al.|1991 HyperLeda, 
[Paturel et al.||2003l GHASP VH; and Muhoz-Mateos et al. 
in prep). We then derived position-velocity diagrams (PV 
diagrams) along the kinematic major axis of the galaxies in 
order to see the rotation in the spatial direction (presented 
as supplementary material in the Appendix G. 

Then, we proceeded with ROTCUR as usual: (1) All pa¬ 
rameters were left free. (2) To obtain the position of the 
galaxy centre, the PA, inclination, and systemic velocities 
were fixed, and the centre values were free. In the majority 
of the cases, the fits were unsatisfactory due to the patchi¬ 
ness of the data. In the galaxies whose centre is represented 
by a bright point source (nucleus), we adopted that as the 
dynamical centre. In the cases where the fits were unsatis¬ 
factory and the centre is not presented with a point source, 
we adopted the value from NED. In any case, the centre po¬ 
sitions differ from those given in NED by less than 0V8, less 
than the angular resolution and less than the accuracy of 
the astrometry in our data cubes. (3) The centre position, 
PA, and inclination were fixed in order to fit the systemic 
velocity. (4) The PA and inclination were then left free with 
all other parameters fixed. (5) Einally, the rotation curve 
was obtained by fixing all parameters except the rotation 
velocity. 

In all the steps, we avoided the sector under 30° from 
the minor axis and a |cos(^)| weight was applied to the 
line-of-sight velocities during the fittings. The objective was 


to minimize the errors caused by the points close to the 
minor axis, where important projections effects occur and 
where the circular velocity term is fitted with difficulty 
[cos{0) ^0]. To determine (fix) one parameter, we forced 
the mean of the values of the rings to satisfy the following 
conditions, based on our experience and linked to the res¬ 
olution of our data: (i) The number of points in each ring 
should be greater than a sufficient number, usually 20-50 
(20 is the minimum number of points in a region with the 
size of our angular resolution, which is arcsec). (ii) The 
difference between the average of all the values from all the 
rings and the value in the given ring should not be more 
than 15 km s“^ for Vsys and 15° for the PA and i (15 km 
s“^ is our velocity resolution, twice the velocity sampling), 
(iii) The uncertainty coming from the free values in the fit 
(i.e., Vsys and the rotational velocity in step 2) should not be 
more than 10 km s“^ for the velocities and 10° for the PA 
and i. Usually, these conditions coincided to pick out those 
cases where the fit was not satisfactory, and the correspond¬ 
ing values were not taken into account in the average. All 
the rotation curves have been computed taking into account 
both approaching and receding sides, and the error bars are 
the dispersion of the rotation velocities computed for all the 
pixels in each elliptical ring. 

To check that the fixed values represent the true rota¬ 
tion, we looked at the residual map (see the following sec¬ 
tion) to search for systematic errors. Eor example, if the 
residual map tends to have negative values, it is because 
the fixed systemic velocity is too large. Other systematic ef¬ 
fects that appear in the residual map when one parameter is 
badly selected are presented in Eig. 8 of [Warner, Wright 


Baldwin (1973). Also, we superimposed on the PV diagrams 


along the major axis the deprojected rotation curve at the 
corresponding angle, confirming that the derived rotation 
curves are correct (see Appendix G. 

We have created two types of rotation curves: a high- 
resolution curve (with a separation of 1” between points) 
and a low-resolution curve (with a separation of 5”). The 
high resolution curve traces better the behaviour of the curve 
and the small-scale features of the galaxy, although it is 
noisier as there are fewer pixels taken into account in the fit 
for each radial point. 

The parameters (z, PA and Vsys) resulting from the low- 
resolution fit are presented in Table Eor NGG 2748, the 
high-resolution rotation curve is presented in Eig. The 
resulting high-resolution rotation curves for all the sample 
galaxies are presented in the Appendix B. 


5.2 Non-circular motions and residual velocity 
fields 

With the rotation curves we explore the circular rotation of 
the sample galaxies. However, there can be deviations from 
these circular motions provoked by dynamical features of 
the galaxy, such as the influence of the potential of the bar, 
past interactions with a companion or streaming motions 
across the spiral arms. To understand the influence of these 
features on the galaxy kinematics, we want to study these 
deviations from rotation, in other words, the non-circular 
motions. 

Eollowing the technique described in Paper I, the first 
step is to create a velocity model map that reflects the ro- 
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NGC 2748 



Figure 2. Results from the analysis of the Fabry-Perot data cubes of the spiral galaxy NGC 2748.Ho; intensity map. b) Velocity map. 
c) Velocity dispersion map. d) Velocity model map. e) Non-circular motions map. f) Ha narrow-band image from AC AM. In the images, 
North is up and East to the left. All the images have the same scale. Similar panels for all the sample galaxies are presented in App. A. 
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Figure 3. High-resolution rotation curve derived from the 5-cr 
velocity map of NGC 2748. Rotation curves for all the sample 
galaxies are presented in the Appendix B. 


tational velocities, which is done by translating the rotation 
curve into a 2D velocity map. The velfi task in GIPSY was 
used for this, assuming kinematic symmetry and using the 
low-resolution rotation curve and the values of i, PA, xo, yo 
and Vsys- To extract the non-circular motions we subtract the 
velocity model map from the observed velocity field, result¬ 
ing in a residual map that is interpreted as the non-circular 
motions map. This residual map shows the deviations from 
pure rotational velocity. In panels d) and e) of Fig. we have 
shown the velocity model and the non-circular motion maps 
respectively for NGC 2748. For the other sample galaxies, 
similar figures are found in the Appendix A. 


Our procedure consists of studying separately those 
parts of galaxies most directly affected by the spiral den¬ 
sity waves (spiral arms), those affected by the potential of 
the bar (bar region), and those affected by both potentials: 
the start of the spiral arms close to the bar (start-of-arms- 
region, SAR hereafter). We have defined these three regions 
(arms, bars and SAR) using images from N ED (mainly SDSS 
false-colour images by Baillard et al.|20n ), 3.6 ym S^G im¬ 
ages (using the output from the 2D decompositions of the 
3.6 /xm images by Salo et al. (2015), to better constrain the 
bars) and our Ha AC AM images. Thus, we are not biased 


© 2015 RAS, MNRAS 000, [^?? 









































Ha kinematics of S^G spiral galaxies-II. Data description and non-circular motions 9 


Table 4. Results. Column I) Galaxy name. Column II) Systemic velocity from optical observations (HyperLEDA). Column III) Systemic 
velocity derived from our FP data. Columns IV and VI) Disc inclination and PA obtained from the 25.5 mag/arcsec^ isophote from 
the ellipse fitting to the S^C 3.6 fim image (Mufioz-Mateos et al. in prep). Columns V and VII) The same as columns IV and VI but 
derived from our FP data. Column VIII) Maximum velocity derived from the rotation curves. Column IX) Quality flag for the maximum 
velocity: reached in our data (I), probably reached (2), probably not reached (3), not reached (4). Notes: (1) PA of the major axis is 
defined as the angle, taken in anti-clockwise direction between the north direction on the sky and the kinematical major axis of the 
galaxy, defined from 0°to 180°. (2) For NCC 1073, 180° have been added to PAs.e^m for a direct comparison with PAna- (3) The AGN 
contribution has been excluded from the measurements. 


Galaxy name 

'^sySjLeda 

(km S“^) 

'dsySjHa 

(km S“^) 

*3.6/im 

n 

iUa 

n 

n 

(°) 

'dmax 
(km S“^) 

dmax 

flag 

NGC 428 

1162.4 

1150 ± 6 

40.3 

45 ± 9 

113.9 

120 ± 4 

127 ± 5 

2 

NGC 691 

2665.0 

2712 ± 5 

41.3 

41 ± 10 

93.7 

91 ± 4 

245 ± 3 

1 

NGC 864 

1558.9 

1550 ± 4 

44.6 

43 ± 5 

22.0 

25 ± 6 

169 ± 8 

1 

NGC 918 

1502.4 

1507 ± 3 

54.6 

57 ± 4 

157.4 

160 ± 2 

153 ± 3 

1 

NGC 1073 

1209.0 

1203 ± 2 

27.3 

29 ± 1 

I 8 O. 5 G) 

165 ± 5 

80 ± 4 

2 

NGC 2500 

479.0 

540 ± 8 

29.9 

41 ± 2 

62.7 

85 ± 5 

111 ±11 

1 

NGC 2541 

530.0 

575 ± 3 

57.0 

57 ± 4 

165.0 

172 ± 3 

96 ± 6 

2 

NGC 2543 

2469.8 

2483 ± 6 

59.2 

61 ± 8 

37.0 

30 ± 3 

210 ± 3 

1 

NGC 2712 

1815.0 

1858 ± 2 

57.9 

58 ± 5 

3.6 

1 ± 3 

176 ± 4 

1 

NGC 2748 

1461.3 

1544 ± 6 

52.9 

74 ± 2 

38.9 

41 ± 2 

150 ± 3 

1 

NGC 2805 

1730.3 

1766 ± 4 

36.0 

36 ± 2 

144.2 

123 ± 3 

106 ± 7 

3 

NGC 3041 

1400.2 

1424 ± 4 

50.7 

50 ± 5 

94.0 

90 ± 5 

173 ± 7 

1 

NGC 3403 

1239.5 

1282 ± 5 

66.9 

66 ± 4 

74.9 

67 ± 4 

167 ± 5 

1 

NGC 3423 

1000.8 

1019 ± 4 

28.4 

28 ± 6 

45.4 

45 ± 5 

177 ± 6 

2 

NGC 3504 

1523.5 

1550 ± 6 

20.9 

39 ± 5 

138.3 

165 ± 8 

151 ± 8 

1 

NGC 4151 

927.7 

1000 ± 5 

48.1 

21 ± 7 

150.6 

20 ± 3 

266 ± 4 

4 

NGC 4324 

1639.1 

1689 ± 5 

63.3 

65 ± 2 

55.0 

57 ± 2 

163 ± 2 

4 

NGC 4389 

712.9 

718 ± 9 

47.3 

45 ± 2 

98.9 

100 ± 3 

149 ± 3 

4 

NGC 4498 

1656.0 

1541 ± 8 

57.0 

58 ± 8 

140.3 

132 ± 7 

120 ± 8 

2 

NGC 4639 

1003.3 

1025 ± 2 

50.2 

39 ± 4 

128.4 

126 ± 4 

222 ± 3 

1 

NGC 5112 

979.1 

1017 ± 5 

49.2 

49 ± 5 

120.4 

123 ± 5 

130 ± 2 

1 

NGC 5334 

1372.0 

1411 ± 5 

41.5 

42 ± 3 

11.9 

10 ± 6 

166 ± 4 

2 

NGC 5678 

1898.5 

1932 ± 5 

56.4 

60 ± 4 

4.1 

3 ± 3 

280 ± 4 

1 

NGC 5740 

1566.1 

1600 ± 3 

56.8 

57 ± 5 

162.5 

162 ± 2 

198 ± 3 

1 

NGC 5921 

1409.4 

1472 ± 3 

33.2 

40 ± 5 

146.7 

150 ± 6 

142 ± 9 

1 

NGC 6070 

2000.1 

2030 ± 4 

63.8 

60 ± 5 

60.0 

57 ± 3 

231 ± 2 

2 

NGC 6207 

848.4 

869 ± 3 

50.5 

57 ± 7 

19.4 

20 ± 3 

136 ± 6 

2 

NGC 6412 

1330.4 

1342 ± 2 

18.4 

20 ± 5 

76.9 

115 ± 5 

175 ± 7 

2 

NGC 7241 

1447.0 

1407 ± 8 

63.5 

65 ± 8 

18.5 

23 ± 4 

142 ± 2 

3 


by a star formation-based classification when defining the 
regions, because there are barred galaxies that show the bar 
in the mid-IR but not in Ha. 


Firstly, based partly on the SDSS images but mainly on 
the IR images, we have defined the region where the stellar 
bar is, also taking into account the bar position angle and 
and bar length (see Sect. 5.3). Secondly, we have identified 
by eye those regions that are between the spiral arms and 
the bar, and which might be influenced by both (SARs). Fi¬ 
nally, we have identified the spiral arms by looking at the 
intensity contours on the SDSS and 3.6 /xm images, drawing 
by eye their extension on the IR images. In the more floc- 
culent spiral arms, it is not easy to identify the extent of 
the spiral arms, so everything except for the bar and SAR 
is assumed to be spiral arms. The three regions have been 
defined astrometrically with irregular shapes, and can be 
identified at all used wavelengths. In Fig. Dl, we present 
the 3.6 fim. S^G images of the barred galaxies with overlaid 
in red and yellow, a schematic representation of the bar re¬ 
gion and the SAR for each barred galaxy. For completeness, 
we have added NGC 5678, a galaxy which may have a bar 


but which has not been classified as either SAB or SB (see 
Sect. 17.2.2). 


To quantify the non-circular motions, we have studied 
the cumulative distribution function (GDF) of the pixel val¬ 
ues in the residual images. We have adopted the value of the 
95% of the distribution of the residual velocities (in absolute 
value) as a representative value of the overall non-circular 
motion within that region. With this method, we avoid in¬ 
troducing possible caveats: if we had used a histogram, the 
bin size would have changed the adopted value for the resid¬ 
ual velocity; and if we had defined regions of certain physical 
scale, the size of the region would matter. Another advan¬ 
tage of the method is that the correlations and trends would 
not change if we choose 50, 70 or 95%, as we represent the 
distribution of all the values. The uncertainties were esti¬ 
mated via a Monte Garlo simulation of the effect of adding 
random noise to the velocity residuals, the former follow¬ 
ing a Gaussian distribution with sigma equal to the spectral 
sampling. In Fig. we present the statistical distribution of 
the residual velocities of the arm region in NGG 2748. The 
adopted value for the residual velocities there is 23.2 zb 1.8 
km s“^; which corresponds to the representative value in 
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Histogram 



CDF 



Figure 4. Statistical distribution of the residual velocities of the 
spiral arm region of NGC 2748. Top) Histogram distribution of 
the residual velocities found in the spiral arm region of the galaxy 
using a bin size of 1 km s“^ (in absolute value). Bottom) Cumu¬ 
lative distribution function of the absolute values of the residual 
velocities for the spiral arm region of NGC 2748. We indicate 
95% value which we use as a representative overall measure of 
the residual velocities. 


the residual map of NGC 2748 (panel e in Fig.|^. The com¬ 
puted values for the non-circular motions in the bar, SAR 
and spiral arms are presented in Table 

For completeness, we have normalised the residual ve¬ 
locities. To do so, we have computed the corresponding cir¬ 
cular velocity at a radius equal to the bar length, to the 
extent of SAR, and to the last measured point (assuming 
that the spiral arms end there), following the universal rota¬ 
tion curve by Persic & Salucci ( |1991 ). This method includes 
several uncertainties, such as assuming that the circular ve¬ 
locity across the complete region is equal to that at the last 
radial point, or that the region is circular with a radius the 
extent of the structure (bar, SAR or arms). Note that the 
residual velocities have been deprojected. 


5.3 Structural parameters 


To understand the nature of the non-circular motions within 
the regions of the galaxies, we study them in conjunction 
with some of the structural parameters that have been de¬ 
rived for the galaxies of our sample. 

The bar strength is quantified by the torque parameter 
Qb, which describes the maximum amplitude of tangential 


forcing normalised by axisymmetric radial force (Combes & 
Sanders 1981| |Buta Block 2001). Qb values have been 


measured from the torque maps derived from t he 3.6 jum 
S^G images (for a description of the method see Salo et al. 


2010 and Laurikainen Sal^|2002 ), and will be published 


m a compilation of bar strengths for the S^G sample (S. 
Diaz-Garcia et ah, in prep). 

The parameter Qb is defined so that Qb decreases with 
the axisymmetric radial force, hence the presence of the 
bulge is implicitly present in Qb. Due to the bulge influ¬ 
ence, Qb is not the best parameter to show the strength of 
the bar on the gas kinematics. However, we will use those 
Qb values as they are the only ones that measure the bar 
strength that have been derived from the S^G images. We 
have also studied the bulge-to-total ratio (B/T). The val¬ 
ues from B/T have been derived from 2D decompositions 
of the 3.6 fim images (Salo et al. 2015| ). We list the Qb 
and B/T values for the galaxies of our sample in Table 
The bar lengths have been measured from the torque maps 
(Diaz-Garcia et ah, in prep.). The bar lengths are depro¬ 
jected units, measured using the method explained in |Salo| 
|et al. ( |201Q ), except when the fits do not have trustwor¬ 
thy quality, in which the bar lengths are measured visually 
(Herrera-Endoqui et ah, in prep.). 

We have also studied the spiral arm class, as a param¬ 
eter that is related to the spiral arm strength. We have 
obtained the arm class classifications carried out by D. 
Elmegreen and presented in the CVRHS, where E means 
flocculent, M denotes multi-armed and G grand design. Note 
that NGC 4324 has not been presented in the table as it does 
not have an arm class classification. The values for the arm 
classification are in Tabled 


5.4 Star formation rates 

The Ho line traces the emission from massive young stars. 
The advantages of this recombination line are that it traces 
very recent SE (timescales ~6-8 Myr), and that it has lower 
sensitivity to dust attenuation than the UV (although not 
negligible, Cardelli, Clayton Mathis||1989 ). Thus, we can 
exploit the Ho imaging and relate the star formation to 
the galaxy’s kinematics. We derive the star formation rates 
(SERs) derived from the Ho luminosity (Lna) following 
[Kennicutt et al. (2009): 


SER(Moyr“^) = 5.5 x 10“^^L(Ho), 


( 1 ) 


where L(Ho) is the luminosity, calculated as 


L(Ho) [erg/s] = 47rD^(3.086 x 10 ^^)^Fhc, 


( 2 ) 


with D the distance to the galaxy in Mpc (Table 1) and 
Ehq: the flux corrected for Galactic absorption, as taken from 
NED and obtained from the Schlally Einkbeiner| ( |2011| 
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Table 5. Non-circular motions as measured in the bar (column II), SAR (column III) and spiral arms (column IV) of the galaxies 
of the sample. We use the 95% value of the cumulative distribution function which we use as a representative overall measure of the 
residual velocities. Column V) Bar strength (Qb) measured from the torque maps derived from the 3.6 fim images (S. Diaz-Garcia 
et al. in prep). Column VI) B/T from the 2D decompositions to the 3.6 p,m images [Salo et ah] ( |20I5] >. Column VII) Arm class (AC) 
classification from the CVRHS: F=flocculent, M=multi-armed, G=grand design. The refers to non-barred galaxies, to those barred 
ones that do not have Ha emission within the bar or to galaxies that do not have an arm classification. (*) These galaxies do not have 
an arm classification in the CVRHS, but as their arm classification in |Elmegreen Elmegreen| ( |I987| was 5, we have classified them as 
multi-armed (M). 


Galaxy 

' name 

'yRES,BAR 
(km s“^) 

'?^RES,SAR 
(km s“^) 

'yRES,ARMS 
(km s“^) 

Qh 



B/T 

AC 

NGC 428 

35.8 


1.4 

23.5 


2.7 

18.8 


2.7 

0.29 


0.03 

0.002 

F 

NGC 691 


- 


- 



19.1 


2.5 

- 



0.170 

M 

NGC 864 

40.1 


1.8 

45.0 


2.4 

17.0 


2.9 

0.47 


0.07 

0.027 

M 

NGC 918 

24.9 


2.3 

26.7 


0.4 

20.2 


2.3 

0.23 


0.02 

0.008 

M 

NGC 

1073 

14.6 


2.4 

13.6 


3.3 

12.8 


3.0 

0.63 


0.08 

0.000 

M 

NGC 

2500 

20.6 


2.4 

16.3 


3.5 

14.5 


3.0 

0.28 


0.03 

0.002 

F 

NGC 

2541 


- 


- 



18.9 


2.2 

- 



0.000 

F 

NGC 

2543 

63.2 


0.3 

43.2 


0.4 

30.8 


1.7 

0.35 


0.08 

0.165 

G 

NGC 

2712 

48.3 


1.5 

23.3 


1.5 

18.5 


2.6 

0.28 


0.05 

0.170 

M 

NGC 

2748 

44.6 


0.4 

21.5 


2.7 

23.2 


1.8 

0.45 


0.03 

0.034 

M 

NGC 

2805 

11.1 


3.3 

13.9 


2.6 

15.2 


2.7 

0.19 


0.01 

0.002 

M 

NGC 

3041 


- 


- 



22.6 


2.3 

- 



0.043 

M 

NGC 

3403 


- 


- 



19.9 


1.9 

- 



0.000 

M 

NGC 

3423 


- 


- 



13.0 


3.1 

- 



0.055 

F 

NGC 

3504 

41.4 


2.0 

11.4 


3.7 

12.6 


3.4 

0.25 


0.06 

0.364 

G 

NGC 

4151 

11.8 


3.1 

10.1 


3.9 

9.6 ± 4.3 

0.09 


0.02 

0.443 

m(*) 

NGC 

4324 


- 


- 



16.2 


3.2 

- 



0.326 

- 

NGC 

4389 

18.6 


2.6 

26.3 


1.5 

21.4 


3.0 

0.52 


0.06 

0.000 

M(*) 

NGC 

4498 

19.7 


2.6 

16.2 


2.8 

23.9 


1.3 

0.46 


0.07 

0.000 

F 

NGC 

4639 


- 


17.7 


2.5 

15.6 


3.0 

0.26 


0.04 

0.112 

M 

NGC 

5112 

14.0 


2.9 

23.4 


2.7 

19.9 


2.1 

0.62 


0.06 

0.000 

M 

NGC 

5334 


- 


18.1 


2.9 

15.3 


2.8 

0.49 


0.08 

0.001 

F 

NGC 

5678 


- 


- 



72.6 


0.1 

- 



0.037 

F 

NGC 

5740 


- 


21.0 


2.4 

21.1 


1.7 

0.16 


0.03 

0.125 

M 

NGC 

5921 


- 


21.9 


1.7 

16.5 


2.6 

0.34 


0.06 

0.111 

M 

NGC 

6070 


- 


- 



20.7 


1.7 

- 



0.045 

M 

NGC 

6207 

18.6 


2.6 

15.5 


3.0 

17.9 


2.3 

0.21 


0.02 

0.000 

F 

NGC 

6412 

12.7 


2.9 

5.2 ± 6.7 

15.0 


2.6 

0.24 


0.02 

0.005 

M 

NGC 

7241 


- 


- 



18.6 


2.3 

- 



0.000 

M 


recalibration of the |Schlegel, Finkbeiner Davis (1998) dust 
map. Equation [^assumes a Kroupa (2001) stellar IMF with 
a mass range of 0.1-100 Mq, an electron temperature of 
Te = 10^ K and electron density Ue = 100 cm“^, whereas 
variations in Te from 5000 to 20000 K would result in a 
variation of the calibration constant (5.5 x 10“^^) of ~15%. 
Variations of Ue = 100 — 10® cm“® would result in variations 


in the calibration constant below 1% (Osterbrock & Ferland 


2006). This calibration also assumes that over timescales >6 


Myr, star formation needs to remain constant, and no infor¬ 
mation about the previous star formation history is given 
(see |Kennicutt Evans||2012 ). 


The uncertainties in the measurements of the SFR come 
from: mainly from the distance, with a typical 20% un¬ 
certainty estimated for the distance value, basic reduction 
processes (i.e. flat-fielding correction, around 2%), the zero- 
point calibration (3%, in agreement with the value of 2% 
typical for photometric nights), uncertainty in the flux mea¬ 
surement (1%) and continuum-subtraction process (around 
11%). Altogether, taking into account the uncertainties re¬ 


lated to the quality of the image (without considering the 
distance) we estimate an uncertainty of 17%. 


To consider only star formation, we have to take into 
account that some of the galaxies of our sample have nu¬ 
clear activity: NGC 918 (no specific classification; Palumbo^] 
Tanzella-Nitti Vettolani||1983 ); NGC 4151 (Seyfert type 


1.5; 

1 . 0 ; 


Khachikian Weedman||1974 ) and NGC 4639 (Seyfert 
Ho, Filippenko Sargent||1997 ), where all nuclear ac¬ 


tivity classifications have been collected in [Veron- Getty <V| 


Veron (2006). The Ha emission in these cases is not all com¬ 


ing from star formation, so the derived SFRs in these galax¬ 
ies should be understood as upper limits. We have masked 
the nuclear emissions in these three cases and excluded them 
from the integration (very small in the cases of NGC 918 and 
the low-luminosity AGN of NGC 4639), so that the global 
Ha luminosities are not biased due to AGNs. NGC 3504 on 
the other hand is a starburst galaxy ( |Balzano]|1983 ) with a 
large amount of star formation in its centre. 


In Fig.|^ we present the SFR as a function of the mor¬ 
phological T-type (from the CVRHS). Excluding the star- 
burst galaxy (which stands out in the plot as a filled green 
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Figure 5. SFR as a function of morphological T-type. Barred 
galaxies are presented with diamonds, whereas non-barred galax¬ 
ies are presented with asterisks. The starburst galaxy NGC 3504 
is identified with a filled green diamond. 


diamond), we find the same trend as widely found in the 
literature (e.g., Kennicutt]|1998 James et al.||20Q4 and ref¬ 
erences therein): early-type galaxies tend to have low SFRs, 
intermediate-type spirals have the highest SFRs and very 
late-type spirals and irregulars have lower SFRs than the 
intermediate-type ones, although higher than those of galax¬ 
ies of the earliest types. Although our sample is small, we 
confirm that the presence of a bar does not affect the global 


SFR (e.g., Phillips 1993 Kennicutt 

J 1994 J 

|Kennicutt 1998 

and references therein; James et al. 

2004 

|Fisher 2006). A 


Kolmogorov-Smirnov (K-S) test on both the barred and the 
non-barred distributions with respect to the SFR results in 
a P-value of 0.094. This P-value indicates the relationship 
between both distributions, and in this case, it means that 
both the barred and the non-barred distributions with re¬ 
spect to the SFR are very similar. 

From our AC AM images, we have measured the total 
SFRs for the galaxies of our sample, and also the SFRs for 
each of the defined regions. In Tablewe present the result¬ 
ing SFRs (total, bar, SAR and arms). 


6 DATA RELEASE 

All the data discussed here are released publicly with 
the publication of this paper: raw but reduced cubes (de¬ 
rot ated, phase corrected and wavelength calibrated cubes), 
continuum-subtracted cubes, 0th, 1st and 2nd order moment 
maps and continuum-subtracted Ha images. The data are 
available in FITS format through the NED and the Centre 
de Donnees Stellaires (CDS). 


7 DISCUSSION 

In this Section we will discuss critically the data and data 
quality, and the results from our analysis in terms of kine- 
matical parameters and non-circular motions. We will not 
analyse here the rotation curves as they will be studied in 


depth in the forthcoming Paper III (Erroz-Ferrer et ah, in 
prep.). 


7.1 Data quality 

We first discuss how our observing, reduction and analysis 
procedures have led to caveats and constraints on the data 
quality. 


7.1.1 Effects from the nature of the Ha data 

One of the scientific goals of this kinematical survey is the 
study of the inner parts of the rotation curves, exploiting 
the high angular resolution provided by GHaFaS data. To 
do so, we have tried not to change the angular resolution 
when reducing the data, as would have been the case when 
using large spatial smoothing kernels or adaptive binning 
methods. As the data are made public, the user can always 
choose to adopt another smoothing kernel depending on sci¬ 
entific interest. 

Normally, Ha data have a typical patchy appearance 
with many blank spaces in between the Hii regions. As we 
have hardly smoothed the data spatially, this effect is rather 
obvious in most of our data. The data set contains galax¬ 
ies with high signal and less patchy appearance (such as 
NGC 2748 or NGC 5678), and also galaxies with low sig¬ 
nal and patchier appearance (e.g., NGC 691 or NGC 918). 
Consequently, the process of deriving rotation curves is dif¬ 
ficult and leads to unsampled radial regions or bad radial 
sampling. For example, the central regions of the galaxies 
NGC 428, NGC 4151, NGC 4324, NGC 4639, NGC 5334, 
NGC 5740 and NGC 5921 are not sampled and can not be 
studied. 


7.1.2 Kinematical parameters versus photometric 
parameters 

Another caveat (though related to the first point) is the 
determination of the kinematical parameters z, PA or Vsys- 
We explained in Sect. |5.l| how we have derived our rotation 
curves using ROTCUR in GIPSY. Due to the nature of our 
data, it is not easy to determine the kinematical parameters, 
and it is more often reliable to use other types of data, with 
full spatial coverage and full spatial sampling albeit at lower 
angular resolution. Thus, in our fits we have used Hi-derived 
parameters, whenever possible, as initial conditions. We also 
use the photometric information provided by the ellipse fit¬ 
ting analysis on the S^G 3.6 /xm images by Munoz-Mateos 
et al. (in prep), although our results do not always match 
the photometric results (see Tablefor a comparison). 

Regarding the PA, the PV diagrams and the position 
of the minor axis in our velocity maps are good tracers of 
the kinematical PA, so the results from ROTCUR describe the 
gas motion better than the photometric results. There is one 
case with a high discrepancy: NGC 4151 with PA3.6=I50.6° 
and PAHa=20°di 3° (the real difference is the same as with 
PAHa=I80°+20°=200°, that is 49.4°). After an inspection 
of the continuum-subtracted cube, and also by looking at 
the velocity map (panel b of Fig. A8 bottom) we have found 
that the PA is closer to 0° , very different from the 150° de- 
rived from the S^G image. Bosma, Ekers & Lequeux ( |l977 ) 
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Table 6. SFRs measured from the AC AM Ho: images, corrected for Galactic absorption (see text). Notes: The refers to non-barred 
galaxies or to those barred ones that do not have Ha emission within the bar. (*) The AGN contribution has been excluded from the 
measurements. 


Galaxy name SFRtot SFRbar SFRsar SFRarms 

(Mo yr-l) (Mo yr"!) (Mq yr"!) {Mq yr"!) 


NGC 428 

0.431 ± 0.172 

0.019 

± 

0.008 

0.078 

± 

0.031 

0.334 

± 

0.134 

NGC 691 

0.628 ± 0.251 


- 



- 


0.628 

± 

0.251 

NGC 864 

0.867 ± 0.347 

0.082 

± 

0.033 

0.035 

± 

0.014 

0.750 

± 

0.300 

NGC 918 

0.475(*) ±0.190 

0.006 

± 

0.002 

0.004 

± 

0.002 

0.465 

± 

0.186 

NGC 

1073 

0.731 ± 0.293 

0.029 

± 

0.011 

0.018 

± 

0.007 

0.685 

± 

0.274 

NGC 

2500 

0.161 ± 0.065 

0.002 

± 

0.001 

0.006 

± 

0.002 

0.154 

± 

0.062 

NGC 

2541 

0.191 ± 0.076 


- 



- 


0.191 

± 

0.076 

NGC 

2543 

0.379 ± 0.152 

0.060 

± 

0.024 

0.022 

± 

0.009 

0.297 

± 

0.119 

NGC 

2712 

0.209 ± 0.084 

0.063 

± 

0.025 

0.021 

± 

0.008 

0.125 

± 

0.050 

NGC 

2748 

0.521 ± 0.208 

0.018 

± 

0.007 

0.158 

± 

0.063 

0.345 

± 

0.138 

NGC 

2805 

1.041 ± 0.417 

0.003 

± 

0.001 

0.017 

± 

0.007 

1.022 

± 

0.409 

NGC 

3041 

0.730 ± 0.292 


- 



- 


0.730 

± 

0.292 

NGC 

3403 

0.234 ± 0.094 


- 



- 


0.234 

± 

0.094 

NGC 

3423 

0.298 ± 0.119 


- 



- 


0.298 

± 

0.119 

NGC 

3504 

0.931 ± 0.373 

0.687 

± 

0.275 

0.148 

± 

0.059 

0.097 

± 

0.039 

NGC 

4151 

0.25l(*) ± 0.100 

0.025 

± 

0.010 

0.009 

± 

0.003 

0.217 

± 

0.087 

NGC 

4324 

0.052 ± 0.021 


- 



- 


0.052 

± 

0.021 

NGC 

4389 

0.061 ± 0.025 

0.035 

± 

0.014 

0.013 

± 

0.005 

0.013 

± 

0.005 

NGC 

4498 

0.112 ± 0.045 

0.011 

± 

0.004 

0.021 

± 

0.008 

0.081 

± 

0.032 

NGC 

4639 

0.170(*) ±0.068 

0.003 

± 

0.001 

0.045 

± 

0.018 

0.122 

± 

0.049 

NGC 

5112 

0.537 ± 0.215 

0.020 

± 

0.008 

0.099 

± 

0.040 

0.418 

± 

0.167 

NGC 

5334 

0.177 ± 0.071 

0.002 

± 

0.001 

0.018 

± 

0.007 

0.157 

± 

0.063 

NGC 

5678 

0.802 ± 0.321 


- 



- 


0.802 

± 

0.321 

NGC 

5740 

0.273 ± 0.109 

0.005 

± 

0.002 

0.062 

± 

0.025 

0.206 

± 

0.083 

NGC 

5921 

0.974 ± 0.390 

0.014 

± 

0.006 

0.107 

± 

0.043 

0.853 

± 

0.341 

NGC 

6070 

0.946 ± 0.379 


- 



- 


0.946 

± 

0.379 

NGC 

6207 

0.405 ± 0.162 

0.007 

± 

0.003 

0.113 

± 

0.045 

0.285 

± 

0.114 

NGC 

6412 

0.484 ± 0.194 

0.030 

± 

0.012 

0.002 

± 

0.001 

0.452 

± 

0.181 

NGC 

7241 

0.109 ± 0.044 


- 



- 


0.109 

± 

0.044 


discussed the first Hi synthesis maps of this galaxy, along 
with a deep photograph by Arp showing the outer spiral. 
They found that the PAkin=19° =b 4°, much closer to our 
PAna than to PAphot- We confirm Bosma’s conclusion: the 
photometric image detects only the oval distribution of the 
galaxy. 


In the same way, there are other galaxies in the sample 
where the mid-IR morphology results in PAphot different 
from PAkin, but our PAkin agree with previous values of 
PAkin found in the literature. These are presented in Tablej^ 
All the galaxies in Tablej^are barred (NGC 2805 is classified 
as SA in the CVRHS, but is claimed to host a bar in the 
2 D structural decomposition of the S^G images, Salo et al. 


2015). Therefore, the differences between the photometric 


and kinematic PA may well be explained with the presence 
of the bar which, depending on its orientation, influences the 
kinematics in the bar region and causes deviations from pure 
circular motions (see Sect. |7.2| . Also, irregular outer spiral 
arms (such in NGG 6412) can influence on the measurements 
of PAphot, and we conclude that PAkin trace the real PA 
better than PAphot- 


As for the inclination, the largest discrepancies occur 
in NGG 4151 (already discussed), and in NGG 2748 and 
NGG 3504, where the inclinations as derived from the kine¬ 
matic data are substantially higher than those derived from 
the surface photometry. This can be due to the non-uniform 


distribution of the Hil regions as compared to the distribu¬ 
tion of the stars and dust as seen in the S^G 3.6 fim images. 
NGG 2748 shows a bulge at faint light levels, hence the dis¬ 
crepancy in inclination. NGG 3504 has a low inclination, 
seen in the S4G image. Also, NGG 2748 and NGG 3504 
host an outer ring (R’ and R’l, respectively), not seen in 
Ra. 


7.1.3 Comparison with literature data 

In the Introduction, we mentioned the Ra FP kinematic 
surveys of |Garrido et al.|2002|(GHASP),[Daigle et al.|2006a 
(SINGS sample) and Ghemin et al. 2006| (VIRGO sam- 
ple). These surveys have used instrumentation similar to 
GHaFaS: Gigale and FaNTOmM (FabryPerot de Nouvelle 
Technologie de I’Observatoire du Mont MeganticQ. These 
were mounted at different telescopes: the 3.6 m Ganada- 
France-Hawaii Telescope (VIRGO and SINGS), the 3.6 m 
European Southern Observatory telescope (VIRGO), the 
1.93 m Observatoire de Haute-Provence telescope (VIRGO, 
SINGS and GHASP), and the 1.60 m Observatoire du Mont 
Megantic telescope (VIRGO and SINGS). All these data are 
seeing-limited, but taken under observing conditions worse 


^ http://www.astro.umontreal.ca/fantomm 
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Table 7. PA from different sources, measured in degrees and defined from North to East. Column I) PAphot from the ellipse fitting 
to the 3.6 fim. images. Column II) PA^in from our Ha FP data. Column III) PA^in values found in the literature. Column IV) 
References. 


Galaxy 

F-^phot 

PAkin,THIS STUDY 

PAkin,lit 

Reference 

NGG 1073 

180.5 

165 ± 4 

164.6 

England, Gottesman & Hunter 

( 1990l 

NGG 2500 

62.7 

85 ± 5 

85 

GHASP VH 

NGG 2805 

144.2 

123 ± 3 

II4 

GHASP VH 

NGG 3504 

138.3 

165 ± 8 

163 

GHASP VH 

NGG 6412 

76.9 

II5 ± 5 

II5 

GHASP VH 


than ours (less than 2” seeing in the very best cases, but 
with seeing values of 8”-12” in the worst cases). In Table 
we show the different quality characteristics of these surveys 
as compared to our GHaFaS data. 

The reduction techniques used in these surveys are 
based on the University of Montreal Improved 3D Fabry- 
Perot Data Reduction Technique^ and are summarised in 
Daigle et al. ( 2006b| . The most important reduction steps 


that differ from ours are spectral and spatial smoothing. 
The three surveys smooth the data spectrally (that is, in 
the wavelength direction), with three-channel or Hanning- 
filtering smoothing. Also, they use the Voronoi adaptive bin¬ 
ning method ( [Cappellari &: Copin||20Q3 ), which consists of 
gathering as many pixels as necessary in order to achieve 
a threshold value of the signal-to-noise ratio for the con¬ 
sidered bin. Hence, the Voronoi method allows maximum 
spatial resolution for regions with strong emission and good 
signal-to-noise values for regions with low emission. 

On the contrary and with the goal of keeping the high¬ 
est spectral and angular resolution possible throughout, we 
have not smoothed the data spectrally and applied a median 
spatial smoothing of 3x3 pixels. These reduction procedures 
have been chosen following our scientific goals of studying 
the inner parts of the rotation curves and studying stream¬ 
ing motions within the different galaxy substructures. 

In Fig.[^ we compare the velocity fields of three galaxies 
in the overlap of our sample with that of GHASP (NGC 864, 
NGC 2500 and NGG 2543) and two galaxies in the overlap 
with the VIRGO sample (NGG 4498 and NGG 4639). We 
compare the velocity fields as obtained from the GHaFaS, 
GHASP and VIRGO data by different reduction methods: 
ours, as explained in Sect. and following the reduction 
techniques used for the GHASP and VIRGO surveys (three- 
channel spectral smoothing and Voronoi adaptive binning, 
see IGhemin et aL]|2006] and [Epinat, Amram Marcelin| 
2008| ). In these figures, it is possible to compare directly the 


advantages and disadvantages of each reduction procedure. 
On the one hand, we see that one drawback of following the 
reduction steps of Sect. is that the spatial coverage is lim¬ 
ited to the location of the Hll regions, resulting in patchy 
maps. But on the other hand, Voronoi smoothing reduces 
the angular resolution often to that typical of Hi data (6”- 
10”). With the results from the adaptive binning, it is very 
difhcult to perform studies inside the structures of the galax¬ 
ies (such as bars or spiral arms), as they cannot be resolved 
any more. This includes the study of the streaming motions 


^ http://www.astro.umontreal.ca/ odaigle/reduction/ 


(this paper), and of the inner parts of the rotation curves 
(Paper HI). As explained in Sect.[^ we publish our raw and 
basic reduced data, and any worker interested in smoothed 
versions of our data sets, either spectrally or Voronoi, can 
simply download the data and perform their own favourite 
recipes. 

In Fig. we present the velocity field of the central part 
of NGG 2543 to show the differences between the GHaFaS 
and GHASP data, both reduced following the techniques 
explained in the current paper. The differences arising from 
the different pixel scale (0/2 of GHaFaS against 0/68 of 
GHASP) and different seeing conditions (1/3 for our obser¬ 
vations against 3/7 for GHASP) can clearly be seen. 

To see the possible impact of the different reduction pro¬ 
cesses, we also compare the rotation curves of the overlap¬ 
ping galaxies. In Fig. we superimpose the rotation curves 
from GHASP and VIRGO papers onto our rotation curves 
for the galaxies in Fig. [Erroz-Ferrer et ^ (2012) give a 
direct comparison between the rotation curve of NGG 864 
from our data and that of GHASP, so this galaxy has not 
been included in Fig. For completeness, we compare the 
rotation curves of the remaining overlapping galaxies in 
Figs. BI-B4. There are many similarities and the rotation 
curves agree reasonably. The possible discrepancies mainly 
come from the different parameters (I’sys, i and PA) used 
for the fits to the velocity fields. Also, the spatial coverage 
resulting from the different reduction techniques is reflected 
in the rotation curves: smoothing kernels such as Voronoi 
help to highlight the emission in the low-signal regions, but 
compromise the spatial resolution. 


7.2 Non-circular motions 

One of the aims of this kinematical study is to analyse the 
influence of the galaxies’ main structural features on their 
kinematics. In other words, we want to study the kinemat¬ 
ical footprints of the components of the galaxies, traced 
for instance by the non-circular motions. Deviations from 
pure circular motion have been widely studied previously 
in the literature. The first studies of streaming motions in¬ 


duced by the spiral density waves were presented by |Bosma| 


(1978), Visser (1980) and Rots et al. 1 

1990 

) using Hi data; by 

Ichikawa et al. (1985), Glemens (1985 

) and 

Gepa et al. (1992) 

using GO data, and by van der Kruit (1976) and Marcelin, 

Boulesteix V Georgelin 

(1985) usin 

g Ha data. They con- 
cted theoretically: vari- 

firm what Roberts (1969) had predl 


ations of 10 — 30 km s“^ in the velocity are found when the 
gas crosses the density wave. Bar-induced non-circular mo¬ 
tions have also been analysed in the literature, and the first 
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Table 8. Characteristics of the Ho; FP kinematical surveys of |Garrido et al.]|2002| (GHASP), |Daigle et al.||2006a| (SINGS sample) and 
[Chemin et al.|2006| (VIRGO sample), shown along with the characteristics of the data used in this paper. Note) The denotes that there 
are no galaxies with that criteria (VIRGO, GHASP and this paper), or that no information about the seeing was presented (SINGS). 
Information before the reduction processes (i.e., before smoothing). 


Parent survey 

VIRGO 

SINGS 

GHASP 

S'^G (this paper) 

Sample size 

30 

28 

203 

29 


Seeing < I'.'S 

11 

- 

- 

29 

2” < Seeing < 4” 

19 

- 

143 

- 

4” < Seeing < 6” 

- 

- 

45 

- 

Seeing > 6'/0 

- 

- 

15 

- 


Angular sampling 

0'/42 - 1"61 

0"42 - I'.'ei 

o"68 - (y.'ge 

072 

Spatial smoothing 

Voronoi S/N=5 
Gaussian 3” x 3” 

Voronoi S/N=5 

Voronoi S/N=5 

Median 3x3 pix 

Spectral sampling^*) 

7-14 km s“^ 

7-14 km s“^ 

~10km s“^ 

~8 km s“^ 


studies were performed by Peterson et al. (1978), Duval & 


Athanassoula (1983) and Pence et al. (1988) using optical 


data; and by Sancisi, Allen & Sullivan (1979) and Gottesman 


et al. (1984) using Hi data. They found that the isovelocity 


contours tend to go parallel to the bar rather than parallel to 
the minor axis. Also, they found some cases where these de¬ 
viations from the circular motion are more prominent when 
the PAbar is at ~ 45° from the kinematic major axis (in 
agreement with the models by, e.g., |van Albada &: Roberts] 
1981). The nature and characterisation of these non-circular 


motions is a topic still under debate, and many studies of the 


streaming motions have been performed (e.g., Knapen et al. 

2000 

Fresneau, Vaughan & Argyle||2005| |Spekkens & Sell- 

wooc 

2007| IGastillo-Morales et al.||2007| |Shetty et al.||2007 

Tam: 

)urro et al.| 

2008| IGarcia-Burillo et al.||2009| |Sellwood| 

& Sanchez|]2010| 

Meidt et al.||2013||Font et al.||2014b|). 

In Sect. |5.2 

we explained how we computed a non- 


circular motions map, the residual velocity field. These resid¬ 
ual velocity fields depend significantly on the derived rota¬ 
tion, that is, on the rotation curve derived from the observed 
velocity field. Consequently, a rigorous analysis of the out¬ 
put of the ROTCUR task in GIPSY has been carried out. We 
want to study both non-circular motions caused by the spi¬ 
ral density waves (hereafter spiral-induced non-circular mo¬ 
tions) and those caused by the gravitational potential of the 
bar (bar-induced non-circular motions). The galaxies in our 
sample span a wide variety of morphological types and fea¬ 
tures, so a statistical approach to the non-circular motions 
is possible, though not straightforward. 


7.2.1 Bar-induced non-circular motions 

We are interested in analysing kinematically the effects of 
secular evolution by investigating the effects of the presence 
of a bar on the rotation of the galaxies. Our kinematic data 
allow us to study the kinematics and the star formation at 
the same time, with the caveat that we are only observing 
the star-forming regions. Our first goal is to understand the 
connection between the presence of a bar and the kinematics 
within the bar region. Our second goal is to study the link 
to star formation within the bar. 

Bars are usually dominated by old population stars, and 


in some cases do not show Ha emission (e.g., in the SB 
galaxies NGC 4639, NGG 5334 or NGG 5921). From the 
12 SAB and 7 SB galaxies in our sample, the bar region 
is clearly defined in Ha in 11 of them. In most of these 
(NGG 864, NGG 1073, NGG 2500, NGG 2748, NGG 2805, 
NGG 4151, NGG 4389, NGG 4498, NGG 5112, NGG 6207 
and NGG 6412) the bar shows Ha emission, and there are 
non-circular motions along the bar (see Table [^. Although 
the bar is not obvious in Ha, using the method shown in 
Fig. El we can in fact measure non-circular motions in three 
further galaxies: NGG 428, NGG 918 and NGG 3504 (Table 
1^. In the other galaxies (NGG 2543, NGG 2712, NGG 5740 
and NGG 5921) the bar itself does not appear in Ha, but 
the surrounding regions (mostly the SAR) show large non¬ 
circular motions. We have also measured the non-circular 
motions in the SAR in all other barred galaxies, and along 
the spiral arms of all galaxies. NGG 4324 is a lenticular 
galaxy, so there we have measured the residual velocities 
along the ring instead of the arms, finding values of 16.2 di 
3.2 km s“^. 

The PV diagrams along the minor axis presented in 
Appendix G are the first clue to see that there are devia¬ 
tions from circular motion. As illustrated for NGG 864 in 
Paper I, in the absence of non-circular motions, the veloc¬ 
ity profile along the minor axis should be completely flat. 
For most of the barred galaxies in our sample, there are de¬ 
viations from the systemic velocity, presented as peaks or 
troughs in the velocity profile along the minor axis. These 
indicate that some gas along the minor axis is not following 
the rotational pattern of the galaxy. As we see in NGG 864 
(Fig. G5), NGG 3504 (Fig. G6) or NGG 5678 (Fig. G7), the 
deviations from the systemic velocity are found beyond the 
extent of the bar. This implies that the non-circular motions 
outside the bar region (i.e., in the SAR) can be a result of 
the influence of both potentials (bar and spiral arms). 

In Fig. |10| we represent the residual velocities in the 
bar region (and normalised) as a function of the difference 
between the major axis of the bar and that of the kinematic 
major axis, APA. We see that the amplitude of the non¬ 
circular motions here does not increase as APA approaches 
45°, contrary to the expectation. 

Fig. (left) shows the residual velocities in the three 
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Figure 7. Comparison between the data obtained with the GHaFaS instrument and data from GHASP (NGC 864, NGC 2500 and 
NGC 2543) and VIRGO (NGG 4498 and NGG 4639). Golumn I) DSS image. Golumn II) AG AM Ha image. Golumn HI) Velocity 
field from GHaFaS observations, obtained following the reduction procedures described in Sect. Column IV) Velocity field from 
GHASP/VIRGO observations, obtained following the reduction procedures described in Sect.[^ Golumn V) Velocity field from GHaFaS 
observations following the reduction processes of GHASP/VIRGO data (three-channel spectral smoothing and Voronoi adaptive binning). 
Golumn VI) Data from GHASP ( [Epinat, Amram &: Marcehn|2Q0^ and VIRGO dChemin et al.|2Q06| > surveys. Note) For NGG 2543, the 
central region has been marked with a black square and explicitly shown in Fig[^ 


regions -bar (Fig. [^), SAR (Fig. [TT| i) and spiral arms 
(Fig. pJ^)- as a function of the bar strength, quantified by 
the torque parameter Qb- Qb is strongly reacting to the 
bulge: a stronger bulge dilutes the tangential force from the 
bar and lowers Qb- It is not necessarily true that the bar is 
weaker, just that motions are more controlled by the spher¬ 
ical potential of the bulge. Therefore, to distinguish the in¬ 
fluence of the bulge in the bar strength, we have also studied 
the non-circular motions as a function of the bulge-to-total 
ratio (B/T) in the middle plots of Fi g. [TT] bar (Fig, pd^), 
SAR (Fig. pJ^ ) and spiral arms (Fig. m- Also, in all the 
plots, we have presented the galaxies with B/T <0.1 with 
blue asterisks, galaxies with 0.1 < B/T < 0.2 with red trian¬ 
gles, and galaxies with B/T > 0.2 (bulge dominated) with 


squares. We list the Qb and B/T values for the galaxies of 
our sample in Table[^ In plots Fig. 11 (a,b,c,d,e,f,g) we only 
represent those cases where a galaxy has a bar (i.e., Qb /B). 
The same plots as in Fig. have been created for the nor¬ 
malised residual velocities, and are presented in Fig. p^ 


Analysing the left hand plots of Figs. 11 and 12 (a,d,g), 
we conclude that Qb does not correlate with the amplitude 
of the non-circular motions. Also, the middle plots of Figs. 

(b,e,f) show no correlation between the amplitude 


11 


and 


12 


of the non-circular motions and the presence of a bulge. 
Although the few cases with a significant bulge {B/T > 0.1), 
show that the residual velocities in the bar decrease as B/T 
increases, more data would be needed to prove whether the 
bulge is constraining the velocities to remain circular. 
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Figure 9. Ha rotation curves for NGC 2500, NGC 2543, NGG 4498 and NGG 4639. Black diamonds correspond to the high-resolution 
rotation curves derived from GHaFaS data (this paper). Red asterisks and blue triangles correspond to approaching and receding rotation 
curves from GHASP data. Pink squares correspond to VIRGO data. 


We have measured the luminosities from the three re¬ 
gions identified previously in the galaxies of the sample (bar, 
SAR and spiral arms), and derived the SFRs in those regions 
using the equations presented in Sect. |5.4| In panels Fig. [H 
and Fig. [T^, we see a slight tendency that the higher non¬ 
circular motions correspond to star forming bars, or to bars 
that show more Ha emission (note that the scale in SFR is 
logarithmic). This is surprising as in the literature, the high 
shear and shocks in the bar region prevent star formation 
in the models (e.g., Athanassoula||1992|) and observations 


(e.g., [Zurita et al.||20Q4| and [Castiho-Morales et al.||2QQ7| ), 
where star formation is inhibited in the bar region and the 
non-circular motions anti-correlate with the Ha luminosity 
there. This would cause a bias in our study, as we might 
not see in Ha those bars with higher non-circular motions. 
Analysing also panels f and i of Figs. Fig. pTjandp^ we see 
that the SFR and the amplitude of the non-circular motions 
do not correlate. 


We previously stated that some galaxies show Ha emis¬ 
sion along their bar, and others do not. Why does that hap¬ 
pen? In the literature, many studies regarding the star for¬ 
mation along bars have been carried out, but with no clear 
answer, (i) As we stated before, the strong shear and parallel 


shoeks limit the formation of stars (e.g.. 

Athanassoula|1992 

Reynaud & Downes 

1998 

Zurita et al. 

2004 

), because the 

giant molecular clouds can be pulled apart ( 

Downes et al. 


1996 Schinnerer et al. 2002). On the other hand, Hll re¬ 


gions have been found under conditions of high shear stress 
and shocks (e.g., Martin Friedh]|1997| Sheth et ^^[2002 


Zurita fc Perez|2(^ ). Martinet & Friedli (1997) found that 

the highest SFRs along the bar correspond to the stronger 
and longer bars in their sample of isolated late-type barred 
galaxies. 

In Fig. we present the resulting SFR values within 
the bar and SAR of our sample galaxies as a function of Qb, 
a parameter which quantihes the strength of the bar and 
consequently that of shocks/shear. From Fig. 


13 


elude that the bar strength does not determine the pres¬ 
ence or quantity of Ha emission along the bar. In our sam¬ 
ple, we see that the strongest bars (Qb,NGCio73=0.63 and 
Qb,NGC5ii2=0.62) show Ha emission, but further study is 
needed to see if this emission is located in the regions of 
high shocks/shear, (ii) |Garcia-Barreto et al. (1996) found 
that the barred galaxies with Hubble types SBa or earlier in 
their sample did not show Ha emission, probably due to the 
low gas content available for inhow (which is not the case, for 
example, for NGC 4151) so the morphologieal type is not in 
principle the only reason, (iii) Martin & Friedli (1997) also 
found that the two highest SFRs correspond to the galax¬ 
ies which present the strongest signs of reeent interaetion 
or merging (NGC 4731 and NGC 7479). (iv) Sheth et al. 
(2000) came to the conclusion that the stars in NGC 5383 
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Bar strength (Qj 5/7’ log(SFR region) 


Figure 11. Measurements of the non-circular motions. The top plots represent the non-circular motions (residual velocites from the 95 % 
of the CDF) in the bar region, whereas the middle and bottom plots represent the residual velocities in the SAR and arms respectively. 
Left) Residual velocities as a function of the bar strength (represented by the parameter Qb from Diaz-Garcia et al. in prep). The 
galaxies without a bar (i.e., Qb=0) have not been represented. Middle) The same residual velocities as a function of B/T. Right) The 
same residual velocities as a function of the SFR in each of the three regions. Note: In all the plots, we have represented galaxies with 
B/T <0.1 with blue asterisks, galaxies with 0.1 < B/T <0.2 with red triangles, and galaxies with B/T >0.2 (bulge dominated) with 
squares. In the top right plot, we have highlighted with a filled green square the starburst nucleus (NGC 3504). 


may have been formed in the spurs before the gas encoun¬ 
ters the dust lane, and travel ballistically through the shock 
at the dust lane, ionizing the regions located at the lead¬ 
ing side of the dust lane, (v) Verley et al. (2007) suggested 
that the presence of star formation in the central region and 
along bars is a consequence of an evolutive sequenee, that 
goes from galaxies with Ha in the bar to galaxies with Ha 
emission in the circumnuclear region but not within the bar, 
ending in galaxies without any emission either along the bar 
or in a central knot. 


study of the physical conditions in the bar related to the for¬ 
mation of stars at the scales of the molecular clouds. Then, 
it may be possible to understand if a recent interaction or 
merger is causing the star formation of the bar, if the galaxy 
has already used up the available gas to be transformed into 
stars, or if the location of the shocks and shear correspond 
to regions of diminished star formation. 


7.2.2 Spiral-induced non-circular motions 


We conclude that there is no simple answer to the ques¬ 
tion why some bars have current star formation and others 
don’t. The answer is probably a combination of the previ¬ 
ous features. With our sample of galaxies, we find that nei¬ 
ther the Hubble type nor the bar strength determine the 
presence of Ha emission within the bar region, but larger 
samples are needed to establish well-based statistical results 
on the SF within bars. Furthermore, deeper and higher- 
resolution kinematic data (e.g., from the Atacama Large 
Millimeter/submillimeter Array, ALMA) would enable the 


In the bar region, the mechanisms triggering star formation 
may be different from those in spiral arms, due to the dif¬ 
ferent dynamics and shock conditions. In Fig. |14| we show 
star formation within the spiral arms and the amplitude 
of the non-circular motions there as a function of the arm 
class, a parameter which might plausibly be related to the 
“strength” of the spiral arms. We do not find any trend or 
correlation. Elmegreen & Elmegreen (1986) found that the 
SER per unit area does not depend on the arm class, and we 
confirm this. Therefore, we agree with other studies which 
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Figure 12. As Fig. |ll| but now for the residual velocities normalised by the circular velocity at the end of the bar region (top plots), 
SAR (middle plots) and spiral arms (bottom plots). These circular velocities have been estimated using the universal rotation curve from 
[Persic Sz SalnccT] ( |1991| >. 


confirm that the SFR does not depend on the strength of 


the arms (Dobbs & Pringle 2009 Foyle et al. 2010), and 
disagree with those that confirm the opposite (e.g., Seigar 
& James 2002 Clarke & Gittins 2006). The presence and 


magnitude of streaming motions in the arms seem to be a 
local phenomenon, unrelated to the SFR or arm class. 

The high amplitude of the non-circular motions found in 
NGC 5678 (with a representative value of ~70 km s“^ above 
the average 20-30 km s~^ for spiral-induced non-circular mo¬ 
tions) leads us to consider that a spiral cannot be the single 
cause, and that this galaxy could be barred. This galaxy was 


classified as SAB in de Vaucouleurs et al. (1991) and argued 
to be barred by jGanda et al. (2007), but the CVRHS re¬ 
classified it as SA from an S^G mid-IR image and therefore 
does not have a value for Qb- The PV diagram along the 
minor axis (Fig. C7) shows the deviations from the circu¬ 
lar motions, very similar to those created by the potential 
of the bar in NGC 864 (Paper 1). We conclude that the 
non-circular motions in NGC 5678 are an evidence support¬ 
ing the claim that NGC 5678 is a barred galaxy, and the 
high amplitude of those non-circular motions could be bar- 
induced rather than only spiral-induced. Another interpre¬ 
tation could be that there was a recent minor merger. The 
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galaxy is asymmetric in the integrated Ha map, contrary to 
the expectations of a symmetric bar. 


8 CONCLUSIONS 

In this paper we present the data from the kinematical study 
of S^G galaxies which started in Paper I, and reach the 
following conclusions: 

(i) We have completed the observations of our kinemat¬ 
ical study of 29 S^G spiral galaxies of all morphological 
types using the GHaFaS FP instrument (FOV of 3.4 x 3.4 
arcmin). The data have seeing limited angular resolution 
(typical values between 0.6 and 1.4 arcsec) sampled with 
0'.'2 per pixel and a spectral sampling ~8 km s“^. These FP 
data have been observed together with Ha flux-calibrated 
images. 

(ii) To reach our scientific objectives, we have followed 
specific data reduction and analysis procedures that guar¬ 
antee high angular resolution (~1”) kinematic cubes and 
data products. 

(iii) The images, data described and rotation curves in 




































































































































20 


S. Erroz-Ferrer et al. 



0 20 40 60 80 

'phot. ( ) 



PAp,.,. (“) 


Figure 6. Comparison between the parameters inclination and 
position angle derived from the S'^G 3.6 /rm images (i.e. pho¬ 
tometric) and derived with the tilted-ring method to our Ha 
FP data (i.e. kinematical). Left) Comparison of the inclinations 
derived photometrically and kinematically. Right) Comparison 
between the position angles derived photometrically and kine¬ 
matically. Note: Due to our definition of the position angles 
(0° ^ PA < 180°), the difference between PA=0° and PA=179° 
is not 179°, but 1°. 


this paper are publicly released through the NED and the 
Centre de Donnees Stellaires (CDS). 

(iv) We flux-calibrate the GHaFaS data cubes following 
the procedures presented in Paper 1. We conclude that the 
flux calibration can not be automated, and no standard cal¬ 
ibration factors can be extracted for any particular GHaFaS 
filter. Instead, flux calibration is specific to each galaxy and 
needs to be performed by comparison with calibrated Ha 
images, such as our AC AM images. 

(v) We have applied the tilted-ring method to our velocity 
maps to extract the rotation curves. Some caveats arising in 
the nature of these FP data need to be taken into account 
(e.g., intrinsic patchiness of the line emission). 

(vi) We have created non-circular motion maps for all the 
galaxies of the sample. We confirm the presence of these non¬ 
circular motions created by the non-axisymmetric potential 
of the bar along the bar region and at the start of the spiral 



Figure 8. Velocity fields of the central region of NGC 2543 re¬ 
duced following the procedures of this paper from GHaFaS (left) 
and GHASP (right) observations. Note the differences due to the 
pixel scale (072 of GHaFaS against 0768 of GHASP) and seeing 
conditions (173 for our observations against 377 for GHASP). 


arms, with a tendency that the more star-forming bars in¬ 
duce higher non-circular motions. We find that Qb does not 
correlate with the amplitude of the non-circular motions. 
However, our data is biased towards bars with recent star 
formation, where strong shocks and shear may take place 
and star formation may be inhibited. 

(vii) Also, we confirm the presence of non-circular mo¬ 
tions created along the spiral arms, but there is no correla¬ 
tion of the amplitude of these non-circular motions with the 
arm class, a parameter that is related to the arm strength. 
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APA (°) 

Figure 10. Top^Amplitude of the non-circular motions (residual 
velocites from the 95 % of the CDF) in the bar region as a function 
of the difference between the major axis of the bar and that of the 
kinematic major axis, APA. Contrary to what has been reported 
in the literature, the highest streaming motions are not found in 
galaxies with APA~ 45°. Bottom) As the top panel but for the 
normalised residual velocities. 
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Figure 14. Top) SFR in the arms per unit area as a function of 
the arm class classification from the CVRHS. Bottom) Amplitude 
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of the non-circular motions correlate with the arm class. Note) 
NGC 4324 has not been represented in the figure as it does not 
have an arm class classification. NGC 5678 has not been repre¬ 
sented either because the amplitude of the non-circular motion 
within its spiral arms presented in Table is understood as being 
caused by the bar (see Sect. |7.2'!^ . 
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